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Abstract. 

We analyze the Levy processes produced by means of two interconnected classes 
of non stable, infinitely divisible distribution: the Variance Gamma and the Student 
laws. While the Variance Gamma family is closed under convolution, the Student 
one is not: this makes its time evolution more complicated. We prove that ~ at least 
for one particular type of Student processes suggested by recent empirical results, 
and for integral times - the distribution of the process is a mixture of other types of 
Student distributions, randomized by means of a new probability distribution. The 
mixture is such that along the time the asymptotic behavior of the probability density 
functions always coincide with that of the generating Student law. We put forward 
the conjecture that this can be a general feature of the Student processes. We finally 
analyze the Ornstein-Uhlenbeck process driven by our Levy noises and show a few 
simulation of it. 



AMS classification scheme numbers: 60E07, 60G10, 60G51, 60J75 
1. Introduction 

Since a few years tlie Levy processes enjoy considerable popularity in several different 
fields of research from statistical physics to mathematical finance (Paul and Baschnagel 
1999, Mantegna and Stanley 2001, Barndorff-Nielsen et al 2001 and Cont and Tankov 
2004 are just a few examples of books reviewing the large body of literature on 
this subject). In the former field, however, the interest has been generally confined 
to a-stable processes which are an important particular sub-class of Levy processes 
(Bouchaud and Georges 1990, Metzler and Klafter 2000, Paul and Baschnagel 1999, 
Woyczyhski 2001), while studies about non stable, infinitely divisible Levy processes 
abound mainly in the latter field (see for example Cont and Tankov 2004 and references 
quoted therein). The appeal of the a-stable distributions is justified by the properties 
of scaling and self-similarity displayed by the corresponding processes, but it must 
also be remarked that these distributions show a few features that partly impair their 
usefulness as empirical models. First of all the non gaussian stable laws always have 
infinite variance. This makes them rather suspect as a realistic tool and prompts the 
introduction of truncated stable distributions which, however, are no longer stable. 
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Then the range of the x decay rates of the probabihty density functions can not exceed 
x~^, and this too introduces a particular rigidity in these models. On the other hand 
the more general Levy processes are generated by infinitely divisible laws and do not 
necessarily have these problems, but they can be more difficult to analyze and to 
simulate. Beside the fact that they do not have natural scaling properties, the laws 
of their increments could be explicitly known only at one time scale. In fact their time 
evolution is always given in terms of characteristic functions, but the marginal densities 
may not be calculable. This is a feature, however, that they share with most stable 
processes, since the probability density functions of the non gaussian stable laws are 
explicitly known only in precious few cases. 

The need to go beyond the processes generated by stable distributions stems also 
from other recent advances in the field of the fractional differential equations. The 
evolution equations of the Levy processes can be put in terms of pseudo-differential 
operators whose symbols are just the characteristic exponents of the processes (Jacob 
and Schilling 2001, Cont and Tankov 2004). The most popular form taken by these 
equations is that of the fractional differential equations, and this generalization of 
the diffusion equations can be put in connection with Levy noises with non-Gaussian 
stable distributions (Gorenflo and Mainardi 1998a, 1998b, Metzler and Klafter 2000). 
It has been put in evidence in a few papers (Chechkin et al 2003, 2004), however, 
that in the case of Levy fiights confined by symmetric quartic potentials the stationary 
probability density functions show two unexpected properties: in fact not only they 
are bimodal, but they also have a finite variance, differently from what happens to 
the non-Gaussian, stable law of the system noise. This suggests that, under particular 
dynamical conditions, the stochastic evolutions produced by stable Levy noises end up 
in non stable distributions, and hence hints to a new physical interest beyond the pale 
of the stable laws. 

Some new applications for the Levy, infinitely divisible but not stable processes 
begin also to emerge in other physical domains (Cufaro Petroni et al 2005, 2006, 
Vivoli et al 2006): as we will see in the following the statistical characteristics of some 
recent model of the collective motion in the charged particle accelerator beams seem 
to point exactly in the direction of some kind of Student infinitely divisible process. 
At the present stage of our inquiry the proposed model for the particle beams is only 
phenomenological and it lacks a complete, underlying, physical mechanism producing 
the noise. This however brings to the fore the problem of the dynamical description 
of complex systems. The infinitely divisible Levy processes with a jump component 
are indeed interesting also in the light of the connection established between Markov 
processes and quantum phenomena by the stochastic mechanics. This latter is a model 
universally known for its original application to the problem of building a classical 
stochastic model for quantum mechanics (Nelson 1967, 1985, Guerra 1981, Morato 
1982, Guerra and Morato 1983), but in fact it is a very general model which is suitable 
for a large number of stochastic dynamical systems (Albeverio, Blanchard and H0gh- 
Krohn 1983, Paul and Baschnagel 1999, Cufaro Petroni et al 1999, 2000, 2003, 2004). 
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As recently proposed, a stochastic mechanics with jumps driven by a non gaussian 
Levy process could find applications in the physical and technological domain (Cufaro 
Petroni et al 2005, 2006). The presence of jumps could for instance be instrumental in 
building reasonable models for the formation of halos in beams of charged particles in 
accelerators. On the other hand this would not be the first time that Levy processes 
find applications in quantum theory since they have already been used to build models 
for spinning particles (De Angelis and Jona-Lasinio 1982), for relativistic quantum 
mechanics (De Angehs 1990), and in stochastic quantization (Albeverio, Riidiger and 
Wu 2001). 

The standard way to build a stochastic dynamical system is to modify the phase 
space dynamics by adding a Wiener noise B(t) to the momentum equation only, so that 
the usual relations between position and velocity are preserved: 

mdQ{t) ^P{t)dt, dP{t) ^ F{t) dt + P dB{t) . 

In this way we get a derivable, but non Markovian position process Q(t). An example 
of this approach is that of a Brownian motion in a fluid described by an Ornstein- 
Uhlenbeck system of stochastic differential equations. Alternatively we can add a Wiener 
noise W(i) with diffusion coefficient D directly to the position equation: 

dQ{t) = v^+){Q{t),t)dt + VDdW{t) . 

and get a Markovian, but not derivable Q{t). In this way the stochastic system 
is reduced to a single stochastic differential equation since we are obliged to drop 
the second (momentum) equation. The standard example of this reduction is the 
Smoluchowski approximation of the Ornstcin-Uhlenbeck process in the overdamped 
case. As a consequence we will now work only in a configuration, and not in a phase 
space; but this does not prevent us from introducing a dynamics either by generalizing 
the Newton equations (Nelson 1967, 1985, Gucrra 1981), or by means of a stochastic 
variational principle (Gucrra and Morato 1983). From this stochastic dynamics, which 
now notably enjoys a measure of time-reversal invariance, two coupled equations can 
be derived which are equivalent to a Schrodinger equation, prompting the idea of a 
stochastic foundation of quantum mechanics. In fact the stochastic mechanics can be 
used to describe more general stochastic dynamical systems satisfying fairly general 
conditions: it is known since longtime (Morato 1982), for example, that for any given 
diffusion there is a correspondence between diffusion processes and solutions of this kind 
of Schrodinger equations where the Hamiltonians come from suitable vector potentials. 
The usual Schrodinger equation, and hence true quantum mechanics, is recovered 
when the diffusion coefficient coincides with h/2m, namely is connected to the Planck 
constant. However we are interested here not only in a stochastic model of quantum 
mechanics, but also to the general description of complex systems as a particle beams, 
and to this end it would be very interesting - as already remarked - to be able to 
generalize the stochastic mechanical scheme to the case of non Gaussian Levy noises. 
The road to this end, however, is fraught with technical difficulties, so that a better 
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understanding of the possible underlying Levy noises should be considered as a first, 
unavoidable step. 

In this light the aim of this paper is to study a few examples of non stable, infinitely 
divisible processes, and in particular we will focus our attention on the Student processes. 
Since the Student family of laws is infinitely divisible but non closed under convolution 
the process distribution will not be Student at every time. We will show however 
that, at least in particular cases, the process transition law is a mixture of a finite 
number of Student laws, and it is suggested that this could be a general feature of the 
Student processes. On the other hand it can also be seen that for every finite time the 
spatial asymptotic behavior always is the same as that of the Student distribution at the 
characteristic time scale; and this turns out to be exactly the behavior put in evidence 
by Vivoli et al 2006 in the solutions of the complex dynamical system used to study the 
behavior of beams of charged particles in accelerators. 

We will limit our considerations to one dimensional models without going into the 
problem of the dependence structure of a multivariate process (see for example Cont 
and Tankov 2004), and we will not pretend any completeness or generality: our aim 
is rather to present the features of a few selected processes to gain a deeper insight 
into their possible general behavior. The paper is organized as follows: in Section [2] we 
recall a few, well known facts about the Levy processes and in particular the connection 
between the transition function p{x,t\y,s) and a triplet of functions s), B{y,s) 
and W{x\ y, s) characteristic of a Levy process. We also propose a different simplified, 
heuristic procedure to find the explicit form of A, B and W: a procedure not completely 
general, but which works well enough for the rather regular transition functions discussed 
in this paper. In the Section [3] we analyze the behavior of two families of laws (the 
Variance Gamma, and the Student laws) which are particular limit cases of a larger 
class of infinitely divisible laws: that of the Generalized Hyperbolic laws which received 
considerable attention in recent years (Raible 2000, Eberlein and Raible 2000, Eberlein 
2001, Cont and Tankov 2004, and references quoted therein). Our two families are in a 
certain sense conjugate to each other since the roles of their probability density functions 
and characteristic functions are interchanged. Let us remark here that all the laws that 
we take in consideration in this paper are infinitely divisible, but - with a few notable 
exceptions - not stable. We then pass in Section H] to study the Levy process produced 
by the Variance Gamma distributions: since this class is closed under convolution, it will 
be easy enough to find both the characteristic triplet, and the laws of the increments 
for every value of the time interval. Apparent similarities notwithstanding, the case 
of the Student processes discussed in Section [5] is rather different from the previous 
one. In fact the Student family is not even closed under convolution so that we do not 
have explicit expressions for the transition laws at every time scale. As a consequence 
we will restrict our attention to a subclass of Student processes by choosing particular 
(but not trivial) values for the parameters, and we will get results about (a) the spatial 
asymptotic behavior of the transition functions at every time, (b) the explicit form of 
the transition functions at time intervals which are integral multiples of a characteristic 
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time constant, and (c) the form of the Levy triplet of functions. In particular we will 
find that, at discrete times, the process law turns out to be a mixture of a finite number 
of Student distributions by means of a new kind of time dependent discrete probability 
distribution. We finally discuss in Section [6] some pathwise properties of our non stable 
processes by showing also a few simulations of the Ornstein-Uhlenbeck processes driven 
by our Levy noises, and conclude with some remark about the perspectives of future 
research. 



2. Levy processes generated by id laws 

A Levy process X{t) is a stationary, stochastically continuous, independent increment 
Markov process. It is well known that the simplest way to produce its transition laws is to 
start with a type of infinitely divisible (id) distributions (see Gnedenko and Kolmogorov 
1968, Loeve 1978, 1987 and Sato 1999 for a more recent monograph): if we focus our 
attention on centered laws, a type of these generating laws can be given by the family 
of their characteristic functions {chf) ip{au) with a spatial scale parameter a > 0. The 
chf of the transition law of our stationary process in the time interval [s, t] will then be 

^au,t - s) = [^{au)f-'^^^ (1) 

where T is a suitable constant playing the role of a time scale parameter, while the 
transition probability density function (pdf) with initial condition X{s) = y, P-q.o. 
will be recovered by an inverse Fourier transform 

p(x,t\y,s) = — lim / $(au, t - s) e"'^^"^^" dw 

= — lim / [^{au)f~'^^^e-'^''-y^''du (2) 

and - because of stationarity - will only depend on the differences x — y and t — s. 

The parameters a and T play a role in the scale invariance properties of the process. 
When the generating family of id laws is closed under convolution the transition laws 
remain within this same family all along the evolution, and the changes are summarized 
just in a time dependence of some parameter of the pdf. But in the case of stable 
laws there is more. If for instance - as in the Wiener process - the generating type 
of law is the normal, centered 7\/'(0, a) it is well known that the transition law (with 
y = and s = for simplicity) is just A/'(0, a^Jt/T), namely it is always normal, but 
with a time dependent parameter: the variance, changing linearly with the time as Dt, 
where D = a? /T is the diffusion coefficient. This means that the overall behavior of 
the process is ruled only by Z^, and not by a and T separately. As a consequence the 
particular values of a and T, namely the particular units of measurement, are immaterial 
and we have the scale invariance. This gives to the Wiener process its property of self- 
similarity: no matter at what space-time scale (namely irrespectively to the values you 
give to a and T, provided that D = o? /T keeps the same value) you choose to observe 
the process, the trajectories always will look the same. 
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These properties of the Wiener process are shared by all the other Levy processes 
generated by stable - even non normal - laws, but not in general by the processes 
generated by other, non stable id laws. It must be remarked, however, that all the non 
gaussian stable laws do not have a finite variance, and show a rather restricted range of 
possible decays for large x: features that partly impair a realistic use of them in empirical 
situations. On the other hand families of non stable, id laws can still be closed under 
convolution, as it is for instance the case of the compound Poisson laws V{\i a; x) with 
cM (p{au) — e^[^("")~^l, where x(w) is the chf of the jump distribution. This means 
again that the evolution of the transition law of a compound Poisson process can always 
be summarized in the time dependence of the Poisson parameter as V{Xt/T, a; x), but 
with respect to the Wiener case there are important differences: while all the transition 
laws of a Wiener process belong to the same (normal) type, Poisson transition laws with 
different parameters do not. The normal laws are indeed stable, while the Poisson laws 
are only id, and Poisson laws with different values of A do not belong to the same type. 
Moreover, while a change in the T value can always be compensated by a corresponding 
change of A so that X/T remains the same, the roles of a and T in a compound Poisson 
process, at variance with the Wiener case, remain completely separated and we do not 
have the same kind of self-similarity. 

The less simple case of processes is finally that generated by families of id laws which 
are not even closed under convolution, since in this event the transition distributions 
do not remain within the same family, and the overall evolution can not be summarized 
just in the time dependence of some parameter. As we will sec in the following this 
is far to be an uncommon situation and this paper is mainly devoted to the analysis 
of particular processes of this kind. It must be kept in mind that in this last case 
the role of the scale parameters becomes relevant since a change in their values can no 
longer be compensated by reciprocal changes in other parameters. This means that, 
to a certain extent, a change in these scale constants produces different processes, so 
that for instance we are no longer free to look at the process at different time scales 
by presuming to see the same features. We should remark, on the other hand, that - 
at variance with the stable, non gaussian case - the pdf 's of the id distributions can 
have both a wide range of decay laws for \x\ +00, and a finite variance a^. For these 
Levy processes generated by id laws with finite variance cr^ it is finally easy to see that 
- due to the fact that the process has independent increments - the variance always is 
finite and grows linearly with the time as aH/T: a feature typical of the ordinary (non 
anomalous) diffusions. 

2.1. The decomposition of a Levy process 

The evolution equations of a process driven by a Levy noise can be given either as partial 
integro-differential equations (PIDE) for the transition functions of the process (Loeve 
1978, Gardiner 1997), or as stochastic differential equations (SDE) for its trajectories 
(Applebaum 2004, 0ksendal and Sulem 2005, Protter 2005). In both cases the structure 
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of the evolution is given in terms of some characteristic triplet of functions. For simple 
Levy process of course this triplet will give rise just to its Levy decomposition in a drift, 
a Brownian and a jump term. In this paper we will choose to follow the description 
in terms of PIDE, and it is important to recall how this characteristic triplet is related 
to the transition functions. We will not attempt to give here a complete and rigorous 
survey of the argument, but we will limit ourselves to fix the notation in a rather 
simplified form (see for example Gardiner 1997, but also for a more rigorous approach 
Leandre 1987, Ishikawa 1994, Sato 1999, Barndorff-Nielsen 2000 and Riischendorf and 
Woerner 2002) suitable for the cases that we will analyze. In particular we suppose to 
consider only processes endowed with well behaved pdf 's, so that (apart from an initial 
distribution) the process is completely defined by its transition pdf p{x,t\ y, s). If then 
we define the triplet of functions 

A(y, s) = lim lim -— / (x — y)p(x, s + Atl y, s) dx (3) 

B(y,s) = lim lim -— / (x — y)'^p(x, s + At\ y, s) dx (4) 

= hm^ — , x^y (5) 

it can be seen that the pdf 's of the process satisfy the following (forward) PIDE 
dtp{x,t) = -d,[A{x,t)p{x,t)] + ^dl[B{x,t)p{x,t)] 

+ lim / [W{x\ z, t)p{z, t) - W{z\ X, t)p{x, t)] dz (6) 

the transition pdf being the solution corresponding to the initial condition 
p(x, s+l I/, s) = S{x — y). In the case of stationary processes (as our Levy processes are) 
the transition pdf p{x, t\ y, s) depends on its variables only trough their differences x — y 
and t — s. As a consequence A and B are simply constants, while W{x\ y, s) = W{x — y). 
It is also known that A plays the role of a drift coefficient, while S is a diffusion coefficient 
connected to the Brownian component of the process; finally W{x\y,s), defined only 
for X ^ y, is the density of the Levy measure of the process. The knowledge of the 
characteristic triplet is also instrumental to write down the PIDE (or alternatively the 
SDE) for other processes driven by a Levy noise. 

In order to calculate the characteristic triplet of a Levy process decomposition 
from ([3]), (jlj) and ([5]) we are supposed to explicitly know its transition pdf We will see 
in the following, however, that given the cM 's of an id distribution it is very easy to 
write the chf ([T]) of the process increments, but also that in general it is not a simple task 
to explicitly calculate the transition pdf by the inverse Fourier transform ([2]). We then 
propose here a different procedure to calculate A, B and W directly from the process 
chf which is surely a known quantity for a Levy process, by adding however that at the 
present stage its derivation is only heuristic. To this end let us remark that from ([T]) 
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p(x, s + At\ y,s) = — lim / [ip(au)f'^^ e-'^^'-^^ciu 

so that, by supposing (which is fair for all the cases that we will consider in this paper) 
ip{—oo) = ip{+oc) = 0, we get with an integration by parts 

p{x, s + At\ y, s) 
At 

= / , lim r [^{au)r/^^e-(^-y^du 
2m{x - y)T m^+oo y„^^ ip{au) 

If now we suppose that our functions are regular enough to allow both to exchange the 
two limits for At — > and for M — > +cxd, and to perform the limit for At —>■ under 
the integral, we immediately have 

/M I / \ 
■M ^(au) 

namely with z = x — y 

W{z) = ^ lim r^e—du (7) 

Remark that in ([7]) the limit must be understood in the sense of the distributions, as can 
be easily checked by applying the formula to some well known case (either the Wiener, 
or the Cauchy process). What is most interesting with respect to the Equation (j5]) is 
that now we can calculate W{z) directly from ip{au), without explicitly knowing the 
transition pdf p{x, t\y,s). 

In the same way for A with an integration by parts we have first of all that 



]- / {x ~y)p{x,s + At\y,s)dx 

J\x-y\<e 

[ [lim r [^{au)r/^^^^e-^-(^-y^d^ 



2mT 



dx 



then, if again it is allowed to freely exchange limits and integrals, we have 
1 f 

lim -— / (x — y)p(x, s + At\ y, s) dx 



At-.0 At J\^.y\<, 

a 



2niT 



\x—y\<e 



I'M 

lim 



dx 



a (p'(au) sinue 
lim / — : — ^ du 



inT Af-»+oo J (p(au) u 



and finally 

A{y, s) = A = —— hm lim / ^ ^ ' du (8) 

ZTtT e-»0+ M-»+oo Lp[au) U 

Here it is understood that the two limits (always in the sense of distributions) and the 
integration must be performed in the order indicated since an exchange will produce a 
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trivial ~ and wrong - result. Remark that when (p{au) is an even function (as happens if 
the process increments are symmetrically distributed around zero), then ' (au) / ip{au) 
is an odd function, and hence - since u~^smue is even - we immediately get A = 0. 
This is coherent with the fact that, when the increments are symmetrically distributed, 
then we do not expect to have a drift in the process. 

As for the coefficient B the usual integration by parts gives 

1 



At 



{x — y) p{x, s + At\ y, s) dx 

\x-y\<e 

[^{au)f'/^ {x - y)e-^<^-y^du 



2mT 



\x-y\<t 



lim ^-ry—ji , N 

A/^+oo (p[au) 



dx 



so that by exchanging limits and integrals we get 

lim -— / (x — y)'^p(x, s + At\ y, s) dx 



'\x-y\<€ 



a 



2'KiT 



\x—y\<e 



dx 



a f ^ (o-u) uecosue — smue 
= —= lim / — - — du 

TlT M->+oo J_j^j ip[au) 

and finally our coefficient is 

a f^^ ip'{au) uecosue -sin ue 
B{y, s) = B = — hm lim / du (9) 

Also in this case we see that for our stationary, independent increment process this 
coefficient is a constant independent from the initial coordinates y and s. 

The formulas ([7]), ([H]) and ^ can finally be checked on two well known (stable) 
cases to give the correct characteristic triplets: the Wiener process produced by a normal 
distribution A/"(0, a) with 

A = 0, ^ = ^' W{z)=Q (10) 

and the Cauchy process produced by a Cauchy distribution C(a) with 

A = 0, 5 = 0, iy(^) = _|_ (11) 

Remark as in these two stable cases the elements of the triplet do not depend separately 
on the two (time and space) scale parameters, but only on a combination of them so 
that a change in the time scale can always be compensated by an exchange in the space 
scale (and vice versa): a point giving rise to the scale invariance which in general is not 
reproduced in non stable processes, as discussed at the beginning of this section. 



3. A class of infinitely divisible distributions 



The increment laws of the Levy processes analyzed in this paper are particular (limiting) 
cases of a larger class of distributions, that of the Generalized Hyperbolic (GH) 
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distributions (for their general properties see for example Raible 2000, Eberlein and 
Raiblc 2000, Eberlein 2001, Cont and Tankov 2004, and references quoted therein). 
The GH distributions constitute a five-parameter class of id, absolutely continuous 
laws with the following pdf 's (for x G M) and cM 's 

^ Kx{8^a^ - p^) (5V«' - (/? + 

where A G M, a > 0, /? G {—a, a), S > 0, /i G R, and Ki,{z) are the modified Bessel 
functions (Abramowitz and Stegun 1968). Apparently a and 5 play the role of scale 
parameters, while /5 is a skewness parameter: the pdf is symmetric when (3 = 0. On 
the other hand /i is just a centering parameter: since in this paper our attention will 
be focused on the symmetric, centered laws, we will always choose = and n = 
and we will consider the more restricted (but still large enough) class Q7i{X, a, S) of the 
centered, symmetric GH laws with the following pdf 's and ciif 's 

fcnix) = ..^^"^^ , ^ {aV6^ + x')^'^ K^_i{aV6^ + x^) (12) 



Vgh{u) = ——===— 13 

with A G M, a > 0, and 5 > 0. 

The GH class contains many relevant particular cases, also for limit values of the 
parameters, and its name comes from the fact that it contains as sub-class with A = 1 
that of the Hyperbolic distributions called in this way because the logarithm of their 
pdf is a hyperbola. The GH distributions arc not always endowed with finite momenta: 
this fact depends on the parameter values and must be explicitly assessed for every 
particular case. On the other hand they are all id, and hence they are good starting 
points to build Levy processes. In general, however, they are not stable laws, and in 
fact they are not even closed under convolution: the sum of two GH random variables 
(rv) is not a GH rv. This means not only that the corresponding processes will not be 
self-similar, but also that often it is not easy to find out what the pdf of the process 
looks like even if it is well known at one time. Remark that the GH class is rich enough 
to contain also as a limit case the sub-class of the normal laws A/'(/i, cr). Indeed it can 
be shown that (in distribution) 

lim lim lim ^7i(A, a, 5) = 7\A(0, a) 

5— »+oo A-+— oo a— »0+ 

provided that 5^/|A| — > 2 cr^. In the following we will study the behavior of the processes 
produced by two other particular limit sub-classes that, at variance with the normal 
distributions, are not stable besides a few exceptions. 
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3.1. The Variance Gamma distributions 

The Variance Gamma (VG) laws (Madan and Seneta 1987, Madan and Seneta 1990, 
Madan and Milne 1991, Madan et al 1998) are obtained from Q'H{\, a, 6) in the limit 
for 5^0+. More precisely, since in general 

\V{v){2/z)\ forz/>0, 
K^{z) = K_^{z) ~ <( - logz, for z/ = , z (14) 

ir(|z/|) (2/^)1-^1, forz/<0. 



we have for A > that 



lim idaYKxiSa) = 2''-'ViX) 

and hence the pdf^s of the centered, symmetric VG laws - which constitute the two 
parameters family VQ{X, a) - are 

fvcix) = ^^J"" ^ («|x|)^-^ K,Ma\x\) . (15) 



2^r(A)V2^ 

with A > and a > 0. As for the corresponding cM 's of VQ{X, a) it is readily seen from 
( IT3|) and (fT4|) that they simply reduce to 

It is apparent that a plays the role of a scale parameter, while A classifies the different 
types of VG laws. For A = 1 the pdf 's and cM 's of the VQ{1, a) laws are 

a -a\x\ ' ^ 



2 ' ^ a2 + M2 

so that VQ{\, a) is nothing but the class of the Laplace (double exponential) laws L{q). 
The pdf f lT^ has an elementary form only when A = n + 1 is an integer with n = 0, 1, . . . 
In fact from 

it is easy to see that (with i = n — j) we have 

For A — > Equation ( fT6l) shows also that our VG laws converge in law to a distribution 
degenerate in x = 0. From the asymptotic behavior of the Bessel functions 

K,{z) ~ e-^ |z| ^ +00 (18) 

we immediately see that the asymptotic behavior of the pdf (fT5l) is (alxD'^e""'^', and 
hence the momenta always exist for every A G M. Of course this corresponds to the fact 
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that the cM (1161) is always derivable in m = 0. Since our laws are centered and symmetric 
the odd momenta vanish; as for the even momenta we have by direct calculation 

, 2'^(2A;- 1)!! TiX + k) 

mya{2k) = ^ , A: = 0, 1, 2, . . . 

so that the expectation is always zero, and the variance 
2 2A 

Then it is easy to see that for a given a > the laws VQ{\, \/2\/a) have all the same 
variance cx^ for every value of A, and that for A — >■ +oo they converge in distribution to 
the normal law A/'(0, a). From the cM (1161) we immediately see that the VG distributions 
are id but not stable. It is easy to see, however, that, the sub-families V^(A, a) with a 
fixed value of a are closed under convolution: in fact the sum of two independent rv's 
respectively with laws V^(Ai, a) and VQ{\2, a) is a rv with law VQ{\i + A2 , a), as can 
easily be seen from (fT6|) . This of course does not amount to stability since laws VQ{\, a) 
and V^(A', a) with A 7^ A' are not of the same type. For the sake of simplicity in the 
following we will take a = 1 and we will use the shorthand notation VQ{\) = VQ{X, 1). 

3.2. The Student distributions 

The class of the centered, symmetric Student laws (see Heyde and Leonenko 2005 for 
a recent review) can be considered as conjugate to that of the centered, symmetric VG 
laws in the sense that here the roles of the pdf and chf are interchanged. They are the 
limit for a ^ 0+ of the QTC{X, a, 6) laws with A < 0. By taking the new parameter 
u = — 2A > 0, and recalling that K^{z) = K^i,{z), the pdf and chf of the QH{X,a,6) 
laws become 



fcnix) - 



Vgh{u) 



so that from Equation f|T4|) in the limit for a 0+ we get the pdf and chf of the 
centered, symmetric Student laws T(z/, 6) 



1+1 
2 



fsAx) = — ^ ( ) (19) 



Vst{u) = 2 ^^.^^y^' (20) 

where u > 0, 6 > and B{z,w) is the Beta function (Abramowitz and Stegun 1968). 
Here 6 is the scale parameter, while u classifies the different law types. It is also easy 
to see that for +00 the Student pdf goes to zero as |a;|~'^~^, so that for a given u 
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the moments msrin) exist only ii n < u. When they exist, the odd momenta are zero 
for symmetry, while the even momenta are 

msTi2k) = 6''^i^^^4f^^ fc = 0,1,2,..., 2k 

In particular the expectation exists (and vanishes) for z/ > 1, while the variance exists 
finite for u > 2 and its value is 

4t = — ^ (21) 



As a consequence, for u > 2 and for a given a > 0, the laws T(z/, cry/u — 2) have all the 
same variance a^, and it is easy to show that for u ^ +00 they converge in distribution 
to the normal law Af{0,a). It can be proved that the Student distributions are id 
(this is not trivial at all; see Grosswald 1976a and 1976b, Ismail 1977, Bondesson 1979, 
Pitman and Yor 1981, Bondesson 1992), but that they are not stable, with one notable 
exception: the z/ = 1 case, that of the Cauchy laws T(l, 6) = C{6) which constitute one 
of the better known classes of stable laws with pdf and chf 

OTT 0^ + 

Besides this case - and at variance with the VG - the Student laws are not even closed 
under convolution: this makes the study of the time evolution of a Student process 
a more complicated and interesting business which constitutes a relevant part of this 
paper. For the sake of simplicity in the following we will take 6 = 1 and we will use the 
shorthand notation T(z/) = T(z/, 1). 

4. The VG process 

Levy processes produced by means of VG distributions are simple enough because of 
their closure under convolution. In fact it is easy to see from flTBl) that (taking a = 1 
and T = 1 to simplify the notations) for a VQ{X) law the transition chf of the process 
(with initial time s = and position y = 0) is 

<l>{u,t\\) = [y,yaiu)f = (^-^y (22) 

so that the law of the increment in [0, t] always is a VG law with the parameter evolving 
in time; namely, at every t, we have X{t) ~ VQ{Xt), and hence the corresponding pdf 
is explicitly known at every time and is 

p(x,t\X) = ^ ^ \x\^'-'^K., i(\x\) (23) 

^ ' ^ 2^*r(At)v^ ^ ^ 

Apparently - as in the Poisson case - the laws of the process belong to the VG family 
all along the evolution, but this does not mean that the process is stable since the 
laws of the VG family are not of the same type. In fact, with increasing values of t, 
the distributions of a VG process go throughout all the gamut of the VG family: what 
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changes with A is just the instant when the distribution is simply a bilateral exponential. 
As remarked in the Section [3TT] the pdf fl23l) has an elementary form only for t = ^, I, . . . 
but a great deal of information is available also in the general, non elementary form. In 
particular from (fT4l) and (fTSjl we can study the behavior of the pdf both near the origin 
and in the asymptotic region. For small x we find 



p{x,t\X) 



|^|2At-i^ forO<t<^, 
-log|x|, fort = ^, rj. 
i_ for ^ < + 

27r r(Ai) ' ^"-"^ 2A ^ ' 



namely near the origin the pdf has an integrable singularity for < t < and 
thereafter it takes finite values for t > As for the asymptotic behavior we have 

t\X) ~ e"!"^!, |x| +00 

namely it is a negative exponential times a power. It is apparent then that this 
asymptotic behavior changes with time since the power depends on t; it is however 
always dominated by the exponential so that all the moments exist at every time. 
From Equation (l23ll we can also explicitly calculate the characteristic triplet: 

e-\A 

A = 0, fi = 0, W{z) = X^ (24) 
so that the dimensionless PIDE for the VG process takes the form 
dM^, t) =lun [ A e-l^l P(- + '^t)-P(-^t) dz 

A validation of (1241) comes then from the Levy-Khinchin formula (Loeve 1987) which 
here reads 

\ogipiu) = Mm [ ( e'"" - 1 - w(x) dx (25) 

and which from (fT6|) and fl2^ easily reduces itself to 

r+oo ^—x 

— log(l + M^) = 2 / (cos MX — 1) dx 

Jo X 

a relation which is immediately verified by direct calculation. 



5. The Student process 

Despite their apparent symmetry and analogy with the VG family, the processes 
produced by the Student laws are not so straightforward to analyze (for recent results 
about the Student process see Heyde and Leonenko 2005). The problem is that the 
Student family T(z/, S) is not even closed under convolution, so that it is not easy to 
figure out the general behavior of a Student process with arbitrary u. As a consequence 
we will limit ourselves here to study the particular case of the u = 3 process whose 
features can be fairly understood: this will also give us an insight on the possible 
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general behavior of these Levy processes. It is important to remark, moreover, that 
this particular Student process with = 3 is the present candidate to describe the 
increments in the velocity process for particles in an accelerator beam (Vivoli et al 
2006), and hence its analysis has not a purely academic interest. Let us introduce now 
the following notation for the T(z/, 5) laws and the corresponding processes: for i/ > 
and 5 > 

/W = = (^) " (26) 

SiA^) = (27) 

so that f{x\ v) from now on will be the pdf of the Student law T{u) = T(z/, 1) for 6 = 1. 
In the same way we can introduce the reduced form of the chf 

Ke{\6u\) 



(p{u\ v) = (p{u\ u, 1) 



2fr(|) 

\u\^ Kk{\u\) 



2fr(|) 

Then, by taking T = 1, the transition chf of the Student process for the law T(z/) (with 
initial time s = and position y = 0) is explicitly known and is 

and the corresponding transition pdf is 

p{x,t\ z/) = — / e-^"^<l>(M,t| u)du = — / e-^"^'[(^(M| u)f du 
2vr y_oo 27r 

If we denote as T(z/, (5)-process the Student process such that its law at t = T is exactly 
T(z/, 6) then p{x, t\ u) will be the pdf of a T(z/)-process. In the following we will perform 
our calculations on the reduced, dimensionless quantities only: we can always revert to 
the dimensional variables by means of simple transformations. It is easy to realize from 
the form of $(m, t\ v) that for t ^ O"*" the process approaches a law degenerate in x = 0, 
and that along the evolution of a Student process the marginal p{x,t \ z/) no longer are 
simple Student pdf^s: after all we know that the Student family is neither stable, nor 
closed under convolution. The main problem is then to find an explicit form for the 
transition pdf which by symmetry can be explicitly written as 



-\ t 

1 



r" + 00 

p(x, t\u) = — I cos(mx 



TT 







2 



2fr(|) 



du (28) 



5.1. The Student processes of odd integer index: the z/ = 3 case 

Since the integration in (128!) can not be performed in general we will limit ourselves to 
particular cases. To do that let us remark that the Student chfs have an elementary 
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Table 1. Examples of odd integer order {v = 2n + 1), dimensionless and reduced 
Student laws. 



ly n /(a;|2n+l) (p{u\2n + I) 



2^-1 g-|u| 



1 + 

-71" \ / 

3 1 + e"l"l(l + |M|) 

5 2 A(1 + :e2)-3 e-|«|(i + |y|+ i|y|2) 

7 3 §il + x'r^ e-H(l + |«| + ||«p + ^|«|3) 



form for odd integer values of the parameter u. In fact from Equation (1171) we have for 
u = 2n + 1 with n = 0,1, . . . and with i = n — j 



^ 1 



on+i "P I 1\ 



so that the chf is just an exponential times a polynomial in \u\ (see Table [U for a few 
explicit examples). The first case n = 0, = 1 is just the stable, reduced Cauchy law 
C(l) which produces the well known Cauchy process. We can then look at the explicit 
time evolution of the first non stable case by taking the n = 1, u = 3 law, namely the 
T(3)-process with pdf 



p{x,t\3) = - / cos(Mx)e"*"(l + u)*(iu 

Jo 



By taking then 

1 f + OO 1 z 

Q{a, z) = - e-^"(l + uY~^du = T{a, z) 

n Jq -k z" 

/+00 
e~V-Mw, r(a,0) = r(a) 

where T{a,z) is the incomplete Gamma function (Abramowitz and Stegun 1968), we 
can also write 

p{x, t\3)=^ {Q{t + 1, t + ^x)} = 5? I '""^ } (29) 

This new closed form (12^ of the increment laws of the Student process with u = 3 \s 
now explicitly given for every time t > 0: in the following sections we will try to analyze 
its properties. 

5.2. Asymptotic behavior of the T (3) -process 

Since the Student laws are not closed under convolution we know that p{x, t \ 3) coincides 
with a Student law only for t = 1. A first question is then to check if, that 
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notwithstanding, some important property of the t = 1 distribution is preserved along 
the evolution. In fact we will see in the following that for an arbitrary fixed, finite t > 
the asymptotic behavior of p{x, t\ 3) for large x is always infinitesimal at the same order 
of the original T(3) 

Proposition 5.1 Ifp{x,t\ 3) is the pdf ^23\j of a Student T{'i) -process, then 

t| 3) = + o (IxT"^) , \x\ +00 

for every given t > 0. 

Proof: Let us remember first of all that by repeated integration by parts of the 
incomplete Gamma function we get the following recurrence formula: for a given a > 
and n = 1, 2, . . . 

le^ r(a) 

Rn{a, z) = - — — V[a - n, z) 

-n z"- V\a — n) 

where, from a classical result about this asymptotical expansion (Gradshteyn and Ryzhik 
1980), the remainder Rn{a,i z) is an infinitesimal of order greater than n 

|i?„(a,z)| =0(|z|-"-^), \z\-^+oo 

Then, for a = t + 1 and z = t + ix with an arbitrary but fixed t > 0, we will have in the 
limit |x| +00 

\'R{Ra{t + l,t + ix)}\ < \Rn{t + l,t + ix)\ =0{\X\-^-^) 

Now take n = 4: from the previous expansion and Equation (129|) we have for |x| +oo 
p(x,t|3) =^{Q{t + l,t + ix)} 

1 A r(t + l) ^{{t-ixf+^] 

while from a direct calculation of the real parts we will find that the higher powers 
exactly cancel away from the numerator so that the leading asymptotic term for 



\x\ 



+00 is of the order |x| ^; more precisely we have 



A r(t + l) ^{{t-ixf^^] _ 2tx^ - At\e - 5t + 3)x2 + 2t\2e - 2t + 1) 
^r(t-fc + l) (t2 + a;2)'=+i (t2 + a;2)4 

= f.+o{\xn 

giving finally the statement in our Proposition. □ 

It must be remarked that the previous result is true for an arbitrary finite, fixed time 
t. For diverging t, however, the reduced law of the process approaches a gaussian: let 
X{t) be our T(3)-process with pdf p{x,t\ 3); then we know that 

E[X(t)] = , Var[X(t)] = E [X'^it)] = t 
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so that t ^/'^X{t) is a centered, reduced rv for every t. A simple look at the ciif's will 
then shows that in distribution we have 
X{t) 



A/'(0,1) 



-oo 



since for large values of t and arbitrary fixed u 



Hu/Vt\3)] 



-\u\/Vi 



1 + 



\U\ 



n t 



-OO 



5.3. The T{'i)-process distribution at integer times t 



n 



To understand the time evolution of t|3) we can analyze the forme of this pdf 
for integral values of the time t = = 1, 2, . . . since in this case the distributions have 
explicit elementary expressions. Of course n\ 3) is nothing else than the distribution 
of the sum of n independent T(3) rv's, so that the following proposition can also be 
seen as a new result about the n-th convolution of the law T(3). 

Proposition 5.2 Forn = 1,2, ... we have (within the notations of the present Section) 



p{x, n\3) = ^ f{x\2k + 1, n) 3) 



qn{k\ 3) 



k=0 



-1 



2k + 



j=0 



2k + 1 
j 



-y 

2n J 



where g„(A;|3) is a discrete probability distribution taking (strictly) positive values only 
for k = 1,2, . . . ,n (in particular g„(0| 3) = for every n) and such that 

qn{k\ 3) ^ 1 
n 



k=l 



2k - 1 



Proof: see Appendix A 



□ 



The meaning of the Proposition 15.21 is then that (at least) at integral times t = n = 
1,2, .. . the marginal one-dimensional pdf p{x,n\ 3) of the T(3)-Student process is a 
mixture (convex combination) of Student pdf^s fl2B]) f{x\i',6) with 



• odd integer orders u = 2k + 1 with A; = 0, 1, . . ., 

• integer scaling factors 6 = n, 

• relative weights qn{k \ 3) such that qn{0 \ 3) = 0, so that no Student distribution of 
order smaller than z/ = 3 appears in the mixture. 

In other words they are mixtures of T{2k + l,n) laws. The distributions qn{k\ 3) are 
a new kind of discrete probability laws whose bar diagrams at different times t = n 
are displayed in Figure [H They show how the weight of the higher order Student 
distributions grows with the time, but also that, this notwithstanding, the lowest order 
{u = 3) distribution is always present - albeit with dwindling importance - with a non 
zero weight. We see at once that this new result is coherent with Proposition 15.11 and 
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Figure 1. Mixture weights of the integer time {t — n) components for a Student 
process with u = 3. 



explicitly shows how the asymptotic behavior is kept 



x\ 



all along the time evolution: 



in fact in the mixture representing p{x, n\ 3) the lowest order Student distribution always 
is - albeit with dwindling weight - that with z/ = 3 which asymptotically behaves as 
all the other components in the mixture are instead faster infinitesimals. The 



\x\ 



-4. 



importance of the higher orders, however, grows with the time. This is exactly the 
behavior recently observed in complex dynamical systems used to simulate the behavior 
of intense beams of charged particles in accelerators (Vivoli et al 2006). Due to their 
mutual interactions these particles follow irregular paths, and a statistical analysis shows 
that the distribution of the increments follows an almost gaussian distribution in its 



central part, and a Student T(3, S) distribution on the tails with a 



\x\ 



decay rate. 



This suggests that the beam particles follow a T(3, 6) Levy process which is observed at 
a time scale (At) large when compared to some characteristic time T of the process, but 
finite and fixed so that the increment distribution shows two different regimes (gaussian 
and |x|~^) in the two regions. 

The results presented in Proposition 15.11 and 15.21 that the pdf of a Levy-Student 
process is a suitable finite mixture of other Student pdf's of different types has been 
proved here only in the particular conditions chosen for our demonstration. It suggests 
however a possible generalization: it is fair in fact to put forward the conjecture that 
every Levy-Student process at every time will have a marginal one dimensional pdf 
which is a mixture of other Student pdf's, but not necessarily (as in our particular case) 
of a finite number of odd integer indices Student pdf's. In other words, by keeping 
always the same notation, the pdf p{x, t|z/o) could be a (possibly continuous) mixture of 
Student pdf's f{x\ z/, 5) through a (possibly continuous) distribution gj(z/| z/q). Finally, 
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Figure 2. Plot of the (reduced and dimensionless) Levy densities for a Student with 
V — 2> (solid line) and for a VG (dashed line) process. 

in order to preserve the result of Proposition 15. ![ we could also conjecture that qt{y \ t'o) 
gives probability zero in the mixture to every Student law with v < vq. If some form of 
this conjecture shows up to be true this would determine some new family of randomized 
Student distributions which is closed under convolution. 



5.4- The Levy triplet for a T{3)-process 

We will finally calculate the elements of the Levy triplet for a T(3)-process from the 
formulas (JTj), ([H]) and ([9]). First of all, due to the T(3) law symmetry, we already know 
that A = 0] then we must recall that the chf of the T(3) law is 

<^(m|3) = e-l"l(l + |u|) (30) 

so that by a direct calculation we get an explicit expression of the Levy triplet with 
T = 1 (for details on the derivation see Appendix B ) 

„ , 1 — Izl ( sin Izl ci l^;! — cos Izl si Izl ) , , 
A = 0, B = 0, Wiz) = ^ -^^2 (31) 

where the sine and the cosine integral functions are (Gradshteyn and Ryzhik 1980) 

/•+°°sint , /•+°°cost , 
SI X = — / at , cix = — at 



t ' J, t 

A plot of W{z) is shown in Figure [2] where it is also compared with the analogous 
density fl24l) for a VG process. The behavior of W{z) at the origin and at the infinity is 

2z-^ + o(z-^), z ^ +00. 



W{z) 



In particular remark that near the origin it has the same behavior of the Levy density 
for the Cauchy process in (fTTj) . while it asymptotically behaves exactly as the T(3) 
distribution. We could then also conjecture here that the W(z) function of a generic 
T(i^, (5)-process will always have a z~^ behavior for z — >• 0+, and a l-zj^'^^^ behavior 
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for \z\ +00. From (IHT]) we also see that, always with T = 1, the PIDE for a 
T(3)-process takes the particular form 



dtp{x,t) = lim / W{z) [p{x + z,t) — p{x,t)] dz 

J\z\>e 



'\z\>e 

with W{z) given in (l3T]) . Finally, inspection into the Levy-Khinchin formula (l25l) for 
the cijf (l30ll immediately gives as a byproduct a previously unknown way to calculate 
a non trivial integral: 

2 /""'"°° sin2;ciz — cos^si^ ,^ , , , ,^ , 

-/ (1 - COSM2)rfz = log(l + M ) (32) 



6. Pathwise properties and simulations 

Both the classes of processes analyzed in this paper do not have a Brownian component 
{B = 0) in their Levy decomposition which is in fact reduced to its jumping part and 
has the form (Cont and Tankov 2004, 0ksendal and Sulem 2005) 

X{t) = / zN{t, dz) + lim / zN{t, dz) 

J\z\>l "^^O^ ie<|2|<l 

N{t,U) = N{t,U)-^u{U) 

where [/ is a Borel set ?7 C M, A^(t, U) is the jump measure of the process, namely is the 
number of the (non zero) jumps of size in U occurring in [0, t], and ^{U) = E [A^(l, U)] is 
the Levy measure of the process. In fact N{t, U) is a Poisson process of intensity z/(t/) 
and N{t, U) is the corresponding compensated Poisson process. The function W{x) 
introduced in the previous sections of this paper plays the role of a density for the 
Levy measure in the sense that ^{dx) = TW{x) dx, so that we have all the elements to 
characterize the Levy decompositions of our processes. In particular, due to the nature 
of the singularities of the W{x) functions in x = 0, it is possible to see that both the 
VG and the T(3, 6) processes (as well as the Cauchy process) have infinite activity, 
namely that z/(]R) = +00. In that event we know (Cont and Tankov 2004) that the 
set of jump times of every trajectory is countably infinite and dense in [0,+oo]. This 
property, together with the continuous distributions of the jump sizes, accounts for the 
fact that at first sight the (simulated) samples of both a VG and a T(3, 6) process do 
not look very different from that of a Wiener process, in particular when we compare 
just the free trajectories of these processes. Then to better see the respective pathwise 
characteristics it will be useful to introduce some Levy diffusions, namely the solutions 
of other SDE driven by a Levy process X{t) (Protter 2004, Applebaum 2004, 0ksendal 
and Sulem 2005). If X(t) is a pure jump Levy process, let us consider the Levy diffusions 
Y{t) solution of the SDE 

dY{t) = a{t,Y{t))dt + dX{t) 

dX{t)= I zN{dt,dz) + lim I zN{dt,dz) 

J\z\>l f^0+ie<|2|<l 
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Table 2. The unit variance laws and pdf 's of the increments used in producing the 
samples of Figure [31 



(a) 


(b) 


(c) 


AA(0,1) 

_}_„-x^/2 
V2¥ 


Vg(l,V2) 


mi) 

2 1 



which is nothing else than a deterministic dynamic system y{t) = a{t, y{t)) perturbed by 
a jump noise X{t). The simplest case is that of a linear force a{y) = —ky giving rise to 
non-Gaussian Ornstein-Uhlenbeck (OU) processes (see for example Barndorff-Nielsen 
and Shephard 2001, Cont and Tankov 2004) 

dY{t) = -k Y{t) dt + dX{t) (33) 

The usual, Gaussian OU process, on the other hand, is the solution of a SDE where the 
noise B{t) is completely Brownian with no jump component: 

dY{t) = -k Y{t) dt + dB{t) (34) 

We can compare now the samples of OU-type processes driven either by a Brownian 
noise, or by a pure jump noise as the VG and the T(3, 6) processes. To do that we 
will produce samples of 5 000 steps by using reduced and dimensionless versions of our 
distributions that we will take of unit variance. In particular we will suppose that for 
time intervals At = T the laws of the noise increments are that reproduced in Table El 
Of course the choice of At = T is instrumental because the VG and the Student laws 
have distributions of elementary form only for At = nT with n integer (and particularly 
simple for n = 1), as we have seen in the previous Sections. It is not so easy, on the other 
hand, to produce our pure jump driven trajectories at other time scales, in particular 
for time scales which are fractions of T. At first sight we could think to overcome this 
difficulty by arbitrarily changing the value of T, but we should remember from our 
previous discussion (Section [2]) that our pure jump processes are not scale invariants, so 
that different values of T produce different processes. Examples of simulated samples 
of these processes are produced by discretizing our SDE and are shown in Figure [31 as 
functions of the dimensionless time r = t/T. The parts (a), (b) and (c) show trajectories 
produced by our three different SDE's: while (a) is a typical sample of an OU process 
solution of the SDE (IMl) driven by a normal Brownian motion, the parts (b) and (c) 
display typical trajectories produced by the SDE f[33|) driven by respectively a VG 
noise and a Student noise. The plots are on the same spatial scale and we can see 
the jumping nature of the non gaussian noises from the fact that, while trajectory (a) 
is rather strictly confined inside the region determined by the restoring force —ky, the 
trajectory (b), and above all the trajectory (c) show clearly random spikes going outside 
the confining region. These spikes are produced by the jumps of the driving processes, 
and the fact that the Student (c) spikes are larger than that of the VG (b) case depends 
on the fact that the VG distribution has exponential tails which - albeit longer than 
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Figure 3. Samples of OU-type diffusions (r = t/T): (a) usual OU process driven 
by gaussian Brownian motion; (b) OU-type process driven by a VG Levy noise; (c) 
OU-type process driven by a Student Levy noise; (d) OU-type process with Student 
noise and restoring force of limited range. 



the gaussian tails - are much shorter than the power tails of a Student distribution (see 
also the corresponding asymptotic behavior of the Levy densities W{z) displayed in the 
Sections m and [53j). The size of the spikes can also be put in evidence by cutting the 
restoring force of the SDE's to a finite length, namely by considering the solutions of 

dY{t) = a{Y{t))dt + dX{t) 

( N J for \y\ < g; 

I 0, for \y\ > q. 

In this case the restoring force acts only when the process lies in [— g,!?], while the 
process is completely free outside this region. Hence when the process jumps beyond 
the boundaries m y = ±g it begins to diffuse freely drifting away from the bounding 
region. Occasionally, however, it can also be recaptured by the binding force. All these 
features are represented in the part (d) of Figure [3] which displays the trajectory of a 
Student driven OU-type process with a limited range of the force. To compare it with 
the other two cases we must now look at the different values of q that make an escape 
reasonably likely: while to let an OU gaussian process to escape is necessary to have a 
rather small value of g, evasions are likely in the VG case for larger, and in the Student 
case even for much larger, q values. 



N Cufaro Petroni: Mixtures in non stable Levy processes 



24 



7. Conclusions 



We have studied in this paper a few examples of non stable, infinitely divisible processes, 
and in particular we have explicitly written down their evolution equations and the laws 
of the increments which are the germ of the corresponding markovian evolutions. In 
particular we focused our attention on the Student processes and we presented a new 
explicit form of their transition functions. Since the Student family of laws is infinitely 
divisible, but non closed under convolution the distribution of the corresponding Levy- 
Student process is a Student distribution only at the one particular time. Along the 
evolution, instead, the process distribution is no longer a simple Student distribution. 
We have shown in the previous sections that, this notwithstanding, at least in the 
case of a specific type of Student distribution (with finite variance), and at least in an 
infinite sequence of equidistant time instants the process transition law is a mixture of 
a finite number of Student laws given by means of a new kind of discrete probability 
distribution. This prompts the conjecture that in fact while the Student family is not 
closed under convolution, some family of mixtures of Student distributions can possibly 
be closed. On the other hand, while it is easy to show that for large values of time the 
reduced increment law tends to be normal (as it should be since we are dealing with 
finite variance distributions), we have also emphasized that for a finite (albeit large) 
time the asymptotic behavior always is the same as that of the Student distribution at 
the unit time. This behavior has been put in evidence by Vivoli et al 2006 in their model 
for halo in particle beams, and we have put forward the conjecture that this could also 
be a more general behavior of the Student processes. This last remark is interesting 
also in connection with a possible generalization of the stochastic mechanics that we 
mentioned in the Section [1] and that will be the argument of forthcoming research. 
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Appendix A. Proof of Proposition 15.21 

For a = n + 1 and n = 1,2,... the incomplete Gamma functions have a finite elementary 
expression (Gradshteyn and Ryzhik 1980) so that 

^/ N , , n\ 1 

Q{n + l,z) = —-J{n + l,z) = Y, 



and hence we get 



p(x,n|3) = -^{Q{n + l,n + ix)} = - V , ... "^l 
^ ^ (n - j)! t 
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(n-j)\ (n2 + x2):'+i I 2£ T ^ 

i=0 ^ "'^ ^ ^ 2^=0 ^ ^ 

where it is understood that the second sum is extended to all the integer values of ^ 
such that < 2£ < j + 1, namely: if j is even then C = 0, 1, . . . , |; if j is odd then 
£ = 0, 1, ... , A little manipulation and the use of Equation fl26l) then give 



71 ^ \j J ^ \ 2£ / ^ \m) n \rfi + x 

j=0 ' 2£=0 ^ ' m=0 ^ ' ^ 

" TT 5^ ^ ( 2£ ) S^"^'' 

j=0 ^ 2e=0 ^ ^ m=0 ^ 

^ ^ ' ~ ^ ^^^^ '^^^ - m) + 1, n) 

with j[x\v) defined in ( l26l) . Now by exchanging the order of the last two sums (with the 
previous conventions about the range of the indexes I and m) we have with k = j — m 



1 . , / \ 'I i \ 
p{x, n\3) = -J2[-)^Yl i-^rfi^\ 2^ -m) + l, n) 

^ j=0 V-^/ ^ 2m=0 

..a^.-»4)l('.v)C) 



2£=2m 

i+1 



^,'1 . 1\ 2^-2™(j + -m)! 

X 5 -, J - m+ * ^^-^ ^ 



j=0 ^ 2m=0 

2'-' 27 m!(j-2m + l)! 

^ b/' (2«)' (J - m)!m!(, - 2m + 1)! -"^1 + "> 

^Viy(2n)%J--_,«(j-*)!(2t-j + l)!'^ ' ' 

where it is understood that the second sum extends over all the k values such that 
j — 1 < 2A; < 2j, namely: for odd j we have k = ^-^^ • • • , J, while for even j we have 
k = |, . . . , j. Finally, by exchanging again the sums and by adopting the convention 
that a binomial symbol (?) always is zero whenever the limitation 6 < a is not verified, 
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we have the results of Proposition 15.21 

n 

p{x, n| 3) = ^ f{x\ 2k + 1, n) g„(A;| 3) 

Since the distribution of our Student process is now represented as a hnear combination 
of the Student T{2k + 1, n) pdf 's, p{x, n\ 3) turns out to be a (randomized, Feller 1971) 
mixture, and the coefficient g„(A;| 3) of this combination must satisfy 

n 

qn{k\3)>0, ^g„(A;|3) = l 

fc=0 

with g„(0|3) = for every n, as can be seen by direct calculation. Hence we have 
that qn{k\ 3) is a discrete probability distribution taking non-zero values only for 
k = 1,2, . . . ,n. Finally by remembering that our Student process has zero expectation 
and variance t = n, and taking also into account the Equation ( 12T]) . we can write 

r"+00 

n = 

k=0 

n 2 
A;=0 

SO that we immediately get also the last result in our proposition. 

Appendix B. Derivation of equation ( 1311) 

From ([9]) and (l30l) we have for a T(3)-process that 

1 [^^ —u ue cos ue — sin ue 
B = — lim lim / j — du 

TC e^0+ M^+oo 1 + \u\ 

2 [^°° smue — uecosue 
= — lim / r du 

71" £-^0+ Jn U(l + U) 

= — lim 

71 e^0+ 



/+00 " p + oo 

n| 3) (ix = 3) / x'^ f{x\2k + l,n) dx 

OO 7 n J —OO 





— — (ci e — e si e) sin e + (e ci e + si e) cos e 







where the sine and the cosine integral functions are defined in the text: hence, as for the 
Cauchy and the VG processes, the Brownian part is absent also in this Student process. 
As for the Levy density W{z), from we get 

W(z) = lim / — e"^"" du 

27fiz M^+oo J_j^j 1 + \u\ 

*^ u 



— — - lim / sm.iu\z\) du 

'k\z\ M^+oo Jq 1 + u 

1 + \z\ (cos \z\ si \z\ — sin \z\ ci \z\) 



TTZ'^ 



SO that for our T(3)-process we finally have 0311) . 
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